Fermion Resonances on a Thick Brane with a Piecewise Warp Factor 
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In this paper, we mainly investigate tiie problems of resonances of massive KK fermions on a 
single scalar constructed thick brane with a piecewise warp factor matching smoothly. The distance 
between two boundaries and the other parameters are determined by one free parameter through 
three junction conditions. For the generalized Yukawa coupling ri'^cj)'^'^ with odd k — 1, 3, 5, . . ., the 
mass eigenvalue m, width F, lifetime r, and maximal probability Pmax of fermion resonances are 
obtained. Our numerical calculations show that the brane without internal structure also favors the 
appearance of resonant states for both left- and right-handed fermions. The scalar-fermion coupling 
and the thickness of the brane influence the resonant behaviors of the massive KK fermions. 
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I. INTRODUCTION 

As we know, in order to unify electromagnetism with 
Einstein gravity, Kaluza and Klein first proposed that 
space-times have more than three spatial dimensions []]. 
Later, higher dimensional space-time with large extra di- 
mensions [2-1 G] had been paid more attention; further- 
more, the emphasis shifted to a "brane world" picture 
which opened up a rich and interesting route towards 
solving the long-standing problems like the mass hierar- 
chy and cosmological constant problems in high-energy 
physics. Subsequently, the idea had been applied in the 
standard model [2, 3, 10]. In particular, inspired also 
in [17, 18], it was put forward that at the energy scales 
of standard model the matter fields cannot propagate 
into extra dimensions whilst gravity on the other hand 
can permeate through all dimensions [9, 11, 12]. A new 
wave of research in the field of extra dimensions came 
with the framework of Arkani-Hamcd, Dimopoulos, and 
Dvali who proposed the large extra dimensions model 
(ADD model) [!)] , which lowers the energy scale of quan- 
tum gravity to 1 TeV by localizing the standard model 
fields to a 4-brane in a higher dimensional space-time 
so that the hierarchy problem can be addressed within 
the framework of ADD model. In Ref. [11], the Randall- 
Sundrum (RS) model provides an alternative explanation 
for the gauge hierarchy problem in particle physics based 
on small extra dimensions and a nonfactorizable bulk ge- 
ometry (AdSs). Subsequent progress about brane world 
and extra dimension suggested that the warped metric 
could even provide an alternative to compactification for 
the extra dimensions [12, 13]. 

In RS warped brane world scenarios, the brane is the 
so-called thin brane, which is an infinitely thin object, 
in which the energy density of the brane is a deltalike 
function with respect to a fifth dimension coordinate. So 
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this model is a very idealized brane world model. Re- 
cently, more realistic thick brane models were investi- 
gated in higher dimensional space-time [19-31]. In thick 
brane scenarios, the coupling between gravity and scalars 
should be introduced. By introducing scalar fields in the 
bulk [■'>'!], the modulus can be stabilized in a RS warped 
brane world scenario. In the thick brane scenario, the 
known solutions can be classified into topologically non- 
trivial solutions and trivial ones. A comprehensive re- 
view on the thick brane solutions and related topics is 
given in Ref. [33] and the definition of thick brane was 
given to avoid the problems related to possible different 
understandings of this term. 

Then, it is possible to describe the matter and interac- 
tions localized on the brane by a natural mechanism in 
higher dimensional space-time. One can consider various 
bulk fields on the brane, such as spin scalars, spin 1/2 
fermions, and so forth [34-42]. However, the localization 
of spin 1/2 fermions shows very interesting and impor- 
tant properties. In addition, the localization problem of 
fermions on a kind of topological defect called domain 
wall has been extensively investigated [35-42]. For lo- 
calizing fermions on the thick branes or domain walls, 
the other kind of interactions were introduced in addi- 
tion to gravity. The common interactions include the 
generalized Yukawa coupling between the fermions and 
the background scalar field. Besides, one can consider 
the gauge field [43, 44] localized in our universe. The 
localization problem can also be discussed in the context 
of supcrsymmetry and supergravity [4-5-47]. In differ- 
ent backgrounds, such as vortex background [48-51] and 
general space-time background [52], the localization of 
various fields was investigated. 

The discontinuous Kaluza-Klein (KK) modes with gap 
called bound states and the continuous gapless states 
with rn^ > can be derived [53-56]. Besides, the 
mctastablc KK states of the graviton or fermion with 
finite lifetime can also be obtained in many brane world 
models [22, 35, 57-G8]. Recently, the problem of the 
fermion resonances has been extensively studied in many 
thick brane models [59, 63-66, 68, 69]. In Ref. [59], the 
authors investigated the general properties of localization 
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of fcrmions and scalars in smoothed field-theoretical ver- 
sions of the type II RS brane world model. They reached 
an important conclusion that, if discrete bound states 
are present in the gravity-free case, those become reso- 
nant states in the continuum, while off-resonant modes 
are highly suppressed on the brane. In Ref. [(>'■>], the au- 
thors investigated fermion resonances on the Bloch brane 
constructed with two background scalar fields [70]. The 
problem of the massive resonant KK modes with even- 
parity for both chiralities was studied in depth and their 
appearance is related to branes with internal structure. 
In Ref. [G4], the fermion resonances were studied in the 
context of a de Sitter thick brane; the authors found that 
there exist resonant KK modes in a single brane without 
internal structure. The two conclusions are different be- 
cause in Ref. [(i4] the coupling parameter 77 was varied 
but in Ref. [■''!] ry was fixed. Particularly, in Ref. [GG], the 
authors extended the fermion resonances to multiscalar 
generated thick branes and found the Kaluza-Klein chi- 
ral decompositions of massive fermion resonances are the 
parity-chirality decompositions. In this paper, we con- 
sider the fermion resonances on the thick brane model 
with a piecewise warp factor [G2]. In addition to the con- 
ventional analyses for the mass eigenvalue m, width F, 
lifetime t, maximal probability Pmax, and mass spectra 
of the massive KK fermions, we give an analysis about 
the effects of a free parameter Vq of the brane on the 
fermion resonances. Some new and interesting results 
are obtained and discussed. 

The organization of this paper is as follows. In Sec. II, 
we review the general aspects of the thick brane model 
with a piecewise warp factor. The equations of motion 
and the solution of the background scalar field are ob- 
tained. Besides, taking into account the interaction be- 
tween the fermion and the scalar field by a generalized 
Yukawa coupling, the Schrodingcr-like equations are de- 
rived with the KK decomposition. In Sec. Ill, we provide 
a complete analysis of the fermion resonances on the thick 
brane in detail. Finally, our discussions and conclusions 
are presented in Sec. IV. 



A. The setup of the thick brane with a piecewise 
warp factor 

Let us consider the single scalar field thick brane 
in (4-t-l)-dimensional space-time with a scalar potential 
V{(j)). Specifically, the action describing such a system is 
given by 



S = 



,(1) 



where R is the scalar curvature, g = det (^mat), -W, N = 
0, 1,2,3,4, and k| = SttGs with the five-dimensional 
Newtonian constant of gravitation. We choose the con- 
formally fiat metric and the line-element is assumed as 



guNdx^Ux^ = e^'^^'-\r|^^dx^'dx'' + dz^), (2) 



where e^^ is the warp factor, the metric tensor presents 
signature in the form (— , +, +, +, and z stands for the 
coordinate of the extra space dimension. In the ansatz, 
the warp factor and the scalar field are considered to 
be functions of z only, i.e., A = A{z), = It is 

straightforward to obtain the Ricci tensor and the curva- 
ture scalar: 

R^, = -(3^'2 + A")?7m- , RiA = -4^" , (3) 



R = -A{iA'^ + 2A")e 



(4) 



where the prime denotes derivative with respect to the 
fifth coordinate z. 

The Euler-Lagrange equations can be obtained from 
the action (1) with the metric ansatz (2). After some sim- 
ple variational calculations, the following coupled nonlin- 
ear differential equations are obtained: 



= 3(A'2 - A") , 



(5) 
(6) 

(7) 



In order to obtain a thick brane, we modify the warp 
factor in [G'i] as follows: 



II. REVIEW OF THE THICK BRANE WITH A 
PIECEWISE WARP FACTOR 



In this section, we will review the construction of the 
thick brane with a piecewise warp factor. With the ac- 
tions as the starting points, the field equations for the 
scalar field and the Dirac field are obtained, respectively. 
Especially, inspired by the warp factor in [G2], we con- 
struct a new warp factor and the solutions of the scalar 
field and the scalar potential are obtained. 



A{z) 



|- cos(vn||2|)+2 j 



-\\n[kl{\z\+Pf 



z < 



\z\ > 



(8) 



where Vq, d, and /3 are constants, fco = v/-A/6 [71] and 
the warp factor is a piecewise defined function. For vari- 
ous values of Vq , the corresponding warp factor is plotted 
in Fig. 1(a). We note that the shape of the warp factor 
here is a normal distribution configuration like in many 
other thick brane models. It is worth noting that the con- 
tinuity of the scalar field (j) and its first order derivative 0' 
closely depends on the continuity of the warp factor A, its 
first order derivative A' and its second order derivative 
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A" at the junctions z = ±d/2, which impose the follow- 
ing three junction conditions among the parameters /?, 
/co, Vb, and d: 



-13 



cos 



+ 2 



T^sin(^^ 



2 J 



2 + COS 
2Vo 



l + 2cos(^) 



2 + cos(^) 



(9) 
(10) 

(11) 



From Eqs. (9), (10) and (11), the parameters /3, ko, and 
d can be derived in terms of Vq, namely, 



2arccos(i^) 



, /3 



1.501 45 



fco = 6.270 8VVo . 

(12) 



The above three conditions can ensure the correct Israel 
junction conditions. The similar situation was explained 
in the work [62]. This means that, in order to get a con- 
tinuous scalar field without derivative jumps, the param- 
eters must satisfy the three junction conditions. Thus, 
with the three junction conditions, we obtain a thick 
brane, and the thickness of the brane d is characterized 
by a free parameter Vq. In the thick brane scenario, the 
warp factor is generally expressed in terms of a "smooth 
function" e^^*^^^ of the fifth coordinate z. Here, the sec- 
ond order derivatives of the warp factor are continuous at 
two junctions z = ±d/2 but the warp factor has third or- 
der derivative jumps. In a sense, we can interpret e^"^'^' 
as a smooth curve [72]. So, we call this kind of thick brane 
the thick brane with a piecewise warp factor matching 
smoothly. Notice that the junction points z — ±d/2 are 
located within the first period of the periodic function 
[(cos(-\/T^|2|) + 2)/3] . In this paper, we focus on the 
problem of the fermion resonances of the massive chiral 
KK modes in the background of the thick brane with a 
piecewise warp factor matching smoothly. 



field and the scalar potential are given by: 

sign(z)2iV2g $(d/2) , |z| > | 
$(z) =F((^(z),p)-6n(q/3,^(z),p) 
+An{q,ip{z),p), 



ip{z) — iarcsinh 



1 , fV%z. 
tan ) 



Vim) = 



p= (-53 + 2077)73, g = 2\/7-5, 

243yo [3 cos(2VVb2)+2 cos(v^2)-2] 
[cos(yVbz)-|-2]'' 



(13) 
(14) 

(15) 



\z\>l 



(16) 



where i is the imaginary unit, F [ip(z),p) is the elliptic 
integral of the first kind and Il(q,ip{z),p) is the incom- 
plete elliptic integral. The shape of the scalar field is 
shown in Fig. 1(c). From the figure, we can see that the 
scalar field has a kinklike structure. Within the inter- 
vals \z\ > d/2, the scalar field tends to two constants. It 
is worth noting that we can ensure that the first order 
derivative 0' is zero when \z\ > d/2 [see Fig. 1(d)]. In 
other words, the extradimensional field profile ^(z) has 
no derivative jumps at z = ±d/2. It is easy to check 
that the curvature scalar R{z) and the scalar potential 
V{(f){z)) are continuous ^. 

Here, the energy density of the thick brane is 



too 



.900 - 2 



dC 



,00 



(17) 



where £ is the Lagrangian density of the matter field. 
The distribution of the energy density Tqq along the extra 
dimension z is illustrated in Fig. 2. The energy density 
has no double-peaked structure which is seen as charac- 
terizing a brane with internal structure [28, 63]. 

We now want to argue that, the massless mode of grav- 
ity (i.e., the four-dimensional graviton) can be localized 
on the brane. We consider the gravitational perturbation 
with the metric fiuctuations given by ds^ = e^^'^^\{ri^i, + 
hf^^{x, z))dx'^dx'^ +dz'^). Furthermore, we can decompose 



/i^^(x, z) in the form h^j 



e2^h^^{x)iirn{z), and re- 



quire that hhe a, four-dimensional mass eigenstate mode 

Uh 



where □ = rj'^'^di^d^. The Schrodinger- 
like equation under the transverse-traceless gauge can be 
obtained[12, 19, 22, 71, 73]: 



■ C(^) + VG{z)i:,n{z) = m^i>^{z) , 



(18) 



From the field equations (5), (6), (7) and the warp fac- 
tor (8), together with the junction conditions, the scalar 



^ We thank the referee for pointing out the problem of the discon- 
tinuity of the scalar potential in the first version of the paper. 
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(a) e' 



.2A(z) 



(b) R{z) 





(c) <t>iz) 



id) 4>'{z) 



FIG. 1. The warp factor e^'^'-^^ curvature scalar R{z), scalar field 0(z), and its first derivative 4>'{z). The parameters are set 
to Vo = 0.2 (dashed thicker), 0.5 (solid thicker), and 1 (solid thinner). 




FIG. 2. The energy density Too{z). The parameters are set 
to Vb = 0.2 (dashed thicker), 0.5 (solid thicker), and 1 (solid 
thinner) . 



with the effective potential 
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(19) 



With the warp factor (8), it is easy to check that this is 
a volcano potential like in the RS model. Furthermore, 
there exists a zero mode solution -00 cx; e^"^^^^. It can 
be shown that the graviton zero mode actually is nor- 
malizable, namely, /J*"^ dz|-0oP oc cq + k^^{d/2 + 

with Co = ^ j^^"^ dz[cos{y/VQ\z\) + 2]^ a constant. So 
the massless four-dimensional graviton can be produced 
by a zero mode[12, 19, 20, 22, 24]. 

Here, the observed four-dimensional effective Planck 
scale Mpi can be derived. As we know, the five- 
dimensional Einstein- Hilbert action of gravity can be re- 
duced to the four-dimensional one, namely. 



(20) 



M 



PI 



where is the fivc-dimensional fundamental scale. 
Thus the four-dimensional effective Planck scale can be 
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derived: 



M 



PI 



Mt 



+00 



-00 

f3r„ , 7,-3/ 



(21) 



The four-dimensional effective Planck scale can be ex- 
pressed through the fundamental scale and a finite one- 
dimensional warped volume. 



B. Spin 1/2 fermion fields on tlie brane 

In this subsection, we will review the general aspects 
of the spin half fermions on the brane. Let us consider a 
massless bulk fermion coupled to the background scalar 
by the generalized Yukawa coupling in five-dimensional 
space-time. The action reads 



Sl/2 



£3 



-5 [^T'^'Dm-^ 



77^'F(</))*] , (22) 



where the spin connection a;j\/ in the covariant derivative 



is defined as 



(23) 



(24) 



In five-dimensional space-time, the spinor structure of 
the four component fermions is determined by ~ 
E^^Y^'i with the being the vielbein and {F^, T^) = 
2^MAf ^ Here we use the capital Latin letters M, N, . . . and 
M,N, . . . to label the indices of five-dimensional space- 
time coordinates and the local Lorentz ones, respectively. 
And are the curved gamma matrices and 7*''^ are the 
flat ones. 

With the conformally flat metric (2), the nonvanishing 
components of the spin connection uj]\{ are given by 



(25) 



Then the flve-dimensional Dirac equation can be derived 
from the action (22) 

[-/^d^ + 7' {d, + 2d,A) - 77e^F(0)] * = , (26) 

where 7''9^ is the four-dimensional Dirac operator on the 
brane. 

In what follows, we focus on the fermion field equations 
of the five-dimensional fluctuations, and reduce the Dirac 
fermion held 4* to the four-dimensional effective field on 
the brane. Following the routine in Ref. [40], we carry out 
the general KK decomposition according to the chiralities 
of the fermions, 



^>{x, z) = e ^^Y. [V'Ln(a;)/L„(z) + V'R«(a;)/R„(^)] 



(27) 



where iph„{x) = -j^ipL„{x) and ipRn{x) = j^ipRnix) 
are the left-handed and right-handed components of the 
Dirac fermion fields on the brane, respectively. They 
satisfy the four-dimensional massive Dirac equations in 
the form of 7''(9^-0Lri(a;) = m„'0H„(a;) and 7^9^V'Rri(a;) = 
'nin'ip'Lnix) ■ At the same time, the KK modes fhniz) and 
fB.n{z) of the chiral decomposition about the Dirac fields 
meet the form of the following coupled equations: 

[9,+rye^F(0)] /l„(z) = m„/R„(z) , (28) 
[d, - ve^F{(b)] h.n{z) = -m„h„{z) . (29) 

With the purpose of obtaining the standard four- 
dimensional effective action for the massive chiral 
fermions, the following orthonormality conditions for 
fhn{z) and fnn{z) are needed: 

+00 p+00 

fhmfhndz ~ / fumfRndz = S^nri , 

-00 J ~oo 

+00 

fhrnfRndz = . (30) 

-oo 

Inspecting of Eqs. (28) and (29), we have 



d 
dz 

d_ 
dz 



d_ 

dz 
d_ 

dz 



+ 77e^F(0)j/L = -m2/L,(31) 



We have thus obtained the Schrodinger-like equations for 
the left and right chiral fermions [G3, 64] 

HLfL{z) = m^fi,{z), (33) 
i^R/R(^) = m^fniz) , (34) 

where the corresponding Hamiltonians are defined as 

i?L= (-^ + ^e^^^(0)) [^^+Ve^Pi'l>)) ' (35) 
H.-(-l-Ve-Fi^)) (|-.e-F(^)). (36) 

In view of (31) and (32), it may seem obvious that the 
Schrodinger equations can be rewritten as 



[^d^,+VL{z)]h = m'h, 
[^d^,+Vniz)]fR^m^fR. 



(37a) 
(37b) 



where the effective potentials of the KK modes have the 
explicit expressions 

Vi (z) = [ve^F(cl>)]^ - [ve^F{(j>)] , (38a) 
Vniz) = ive^Fi^)]^ + a,[77e^F(0)] . (38b) 



III. FERMION RESONANCES ON THE THICK 
BRANE 

In this section, we mainly discuss the problem of the 
fermion resonances on the thick brane with a piecewise 
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warp factor. Specifically, we investigate the effects of 
two parameters k and Vq on the resonant states of the 
massive KK modes, respectively. Here /c is a parameter in 
the generalized Yukawa coupling jyvj/^'^^', fc = 1, 3, 5, . . .. 
Our discussions focus on the properties, i.e., the mass 
spectra, the lifetimes, and the number of the massive 
resonant KK modes, related to the two parameters. 

In Rcf. [G2], the parameter x which effectively 
parametrizes the thickness of the domain wall has a very 
interesting impact on the gravitation resonances. In this 
paper, we want to know the impacts of the parameters 
on the behaviors of the fermion resonances. 

In this section, some new natures of the resonant states 
are obtained. We will show that the number of the res- 
onant states increases when k becomes larger, however, 
it decreases as Vq increases. The height of the potential 
function of the resonant fermions with left-handed chi- 
rality is the same as that of the one with right-handed 
chirality. The parameter Vq determines the depth of the 

I 



potential well of the resonant fermions, but it has no in- 
fluence on the height of the potential function. 



A. The fermion zero modes 

As we know, in order to study the localization and 
resonance (quasilocalization ) problems, the coupling be- 
tween the fermion field and the background scalar field 
should be introduced. Here, we choose the generalized 
Yukawa coupling = cf)'' with odd fc = 1, 3, 5, . . .. 

In this subsection, we mainly study the localization of 
the zero modes of left- and right-handed fermions on the 
branc. 

With the given k ~ 1 and brane solution (t>{z) in 
Eq. (13) obtained in the previous section, the potential 
functions can be rewritten in the following specific forms: 



ni 

81 



r]P(j)^{z) + 18(^(z)^Aosin(^/K;I) -f 



Vn{z) = Vi^{z)\„^^r,, 





SCC^ 


m 


+4;-i6 








2)](. 





\z\ < 



z > 



(39a) 



(39b) 



where I ~ 2 + cos{\/Voz). 

In what follows, we analyze the asymptotic behav- 
ior of the potential functions for the left- and right- 
handed fermions. After simple calculation and analysis, 
we find that both the potential functions have a very 
good asymptotic behavior. Specifically, when z tends to 
zero and infinity, their asymptotic values are constants. 
It can be expressed as follows: 

Vi(0) = -Fr(0) = -VMv, (40) 
1/l,r(2 ^ ±w) ^ . (41) 

It can be seen that, for the same set of parameters, 
the potential function of left-handed fermion KK modes 
Vi^{z) has an opposite value compared to the one of right- 
handed fermion KK modes Vr(z) at the coordinate origin 
z = 0. When z tends to infinity, both potential functions 
vanish. The behaviors of the two potential functions are 
illustrated in Fig. 3 for given Vq and various values of rj. 
From Eq. (38) and the continuity of the first derivatives 
of the scalar field and warp factor, we know that the po- 
tentials Vl(z) and Vr(z) are continuous. However, they 
are not smooth because the first order derivative of the 
scalar field is not smooth, this can be seen from Fig. 1(d). 
The insets in Fig. 3 show the enlarged view of the conti- 
nuity of the potentials for clarity. From the left panel of 
Fig. 3, we can see that Vl(z) is a modified volcano-type 



potential like in most of the other contexts. So there is 
no discrete mass spectrum of bound states, and there ex- 
ists no mass gap to separate the massless mode from the 
massive KK modes for Vi^{z). 

In Fig. 3, from the right panel, we see that there is 
a potential barrier around the brane location for Vr(z). 
Nevertheless, for stronger coupling (bigger 77), the poten- 
tial well emerges at the top of the barrier. This means 
that the behavior of the potential function is related to 
the scalar-fermion coupling. In order to clearly examine 
the impact of Vq on the behavior of VFi,(z), we fix k and 
1] in the following. The shapes of Vl{z) and Vr(z) for 
various values of Vq can be depicted easily. As we know, 
for Vr(2:), in general, the potential barrier can not trap a 
massive fermion. However, when Vq becomes smaller, the 
potential well located at the top oiVR{z) is deep enough, 
and a massive fermion with a finite lifetime would appear, 
however, the height of the potential remains unchanged. 
This means that the phenomenon of massive resonant 
states will also occur. 

If we recast the scalar-fermion coupling rj'^cj)'^ {k — 1) 
to 77^(/)'^^' (fc = 3), the configurations of the two poten- 
tial functions Vl{z) and Vb.{z) will change. In this sit- 
uation, Vl(z) becomes a double- well potential, but the 
right-handed one is still a single- well potential. However, 
for the two potential functions, if the other parameters 
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(a) V^{z) 



(b) Vr{2) 



FIG. 3. The potential functions Vi,{z) and Vb.{z) for ^(0) = (j)- The insets show the enlarged view of the junction points for 
clarity. The parameters are set to Vb = 0.5, rj = 0.7 (dashed thicker), 0.5 (solid thicker), and 0.2 (solid thinner). 



are fixed, when the parameter Vq becomes smaller, the 
potential well will become deeper. Furthermore, at the 
center of the brane along the direction of the extra di- 
mension, the values of the potential functions Vi,^^{z) are 
zero, i.e., 14,(0) = Vr(0) = 0. This result is similar to 
that obtained in Refs. [64, 66]. 

As discussed above, Vi,{z) is always a modified vol- 
cano type of potential. It is well known that the poten- 
tial of the left-handed fermions does not provide mass 
gap between the zero mode and KK excitation modes, so 
there is a continuous gapless spectrum of KK excitation 
modes. From Eqs. (28) and (29), the zero mode of the 
left-handed fermions is derived as 

fw{z) cx exp (^-77^ dze^(^")0(z)^ 
r exp (-f dz[cos(^A^|z|) + 2fm) , < f 



dition 



exp 
• exp 



(-f/o'/'df[cos(yT^|z-|) + 2]20(z)) 



Id/2 |2|+/3 



Z > 



When \z\ > d/2, in the first factor, because the integrand 
and the integrating range are finite, the integration result 
must be a constant: 



ci = exp — — 



d/2 



dz[cos{y^\z\) + 2]'^<P{z) . (42) 



In the second factor, we know that (f>{z) = C2 (a constant) 
when z > d/2. So, we obtain 



/Lo(.)ocexp(^-^/^^ 



Id/2'^^\z\+l3^ 
C3(kl+/3)"'^, 



(43) 



where C3 = {d/2 + P) 



riC2/ko 



Then the normalization con- 



+00 



floiz)dz = 1 



(44) 



is equivalent to 



+00 



dz{z + P) '=0' < 00. 



(45) 



d/2 



It is clear that the normalization condition is turned out 
to be 2r]C2/ko > 1. Hence, the zero mode (42) is nor- 
malizable only if rjC2/kQ > 1/2. Besides, although there 
is no explicit expression for /lo(^;), we can numerically 
integrate (42). The plot of the normalizable zero mode 
fhoiz) via numerical integration is given in Fig. 4. From 
Fig. 4, we can better understand the process of the cal- 
culation. On the other hand, the potential VYi{z) around 
the brane location is always positive, and it gradually be- 
comes zero when z — > 00. We know that this type of po- 
tential cannot trap any bound state fermions with right 
chirality and there exists no zero mode of right-handed 
fermions. This result is consistent with the previous well- 
known conclusion that massless fermions must be single 
chirality in the brane world models [3-^), G4]. 

Nevertheless, for a given Vq, the structure of the poten- 
tial Vr is determined by the parameters k and the cou- 
pling constant rj jointly. For a given fc, when rj increases, 
a potential well around the location of the brane would 
emerge and the well would be deeper and deeper. In this 
case, there may exist massive resonant KK modes. In or- 
der to make a potential well appear, we need very large 
77, i.e., the coupling between the scalar and the fermion is 
very strong [64, 66]. In this paper, we would not consider 
this case and our main task is focusing on the impacts of 
parameters k and Vq on the massive resonant states. 
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FIG. 4. The normalizablc formion zero mode fu){z). The 
parameters are set to 77 = 1.0, Vo — 0.2 for the dotted hne 
and Vo = 0.5 for the dashed hne. 



B. The massive KK modes 

As mentioned above, in this subsection, we mainly dis- 
cuss the problem of the fermion resonances on the thick 
brane. Some interesting properties and results will be ob- 
tained. Following the routine in Refs. [()4, (iG], with the 
potential functions (39), we get the numerical solutions 
of the Schrodinger equations (37) with the Numerov algo- 
rithm. Furthermore, the resonant probabilities, lifetimes 
and mass spectra of the resonances will be obtained. 

When the parameter k in scalar-fermion coupling is 
treated as a variable, the maxima of potential functions 
Vh{z) and V^{z) are the same magnitude and the prop- 
erties of resonances are same to the results obtained in 
Refs. [64, GO]. When we treat the parameter Vq as a vari- 
able, the maxima of potential functions Vi^{z) and V^,{z) 
are also the same magnitude but the resonant behaviors 
of the massive chiral fermions with the left- and right- 
handed chiralitics arc different from the case in that k is 
treated as a variable. But the resonant states satisfy the 
Kaluza-Klein parity-chirality decompositions of massive 
fermion resonances obtained in Ref. [{>{>], because the res- 
onant peaks have a one-to-one correspondence between 
left- and right-handed fermions at some resonant mass 
eigenvalue. 



1. Case I: k as a variable 

In this subsection, we mainly study the impacts of 
the parameter fc, which is related to the coupling type 
between the scalar field and the fermion, on the mas- 
sive resonant KK modes. From the context of numeri- 
cal analysis, we know that one should impose the initial 
or boundary conditions for the second order differential 
equations (37). Here, we attach two types of initial con- 
ditions for Eqs. (37). They are 

/(0) = rfo, /'(0)=0, (46) 



and 

/(O) = 0, /'(O) = di . (47) 

Because the potential functions have a Z2 symmetry, the 
solutions of the Schrodinger equations (37) with even par- 
ity and odd parity will be obtained for the above two con- 
ditions, respectively. In this paper, the constants do and 
di arc set to do ^ I, di = 5 for the sake of convenience. 
From the point of view of quantum mechanics, we know 
that the massive KK modes will feel a combined effect of 
the potential barrier and the potential well around the 
location of the brane. Therefore, the massive KK modes 
will show different natures when the depth of the poten- 
tial well changes. 

For the simplest generalized Yukawa coupling type 
F{(j)) = (f>^ . /c = 1,3, when we choose a Vq, some res- 
onant states would be obtained. When the eigenvalues 
deviate from the mass eigenvalues of the resonances, the 
resonant states will become unclear and even disappear. 
The graphics of massive resonant wave functions, poten- 
tial functions, and off-resonant wave functions are shown 
in Fig. 5. From the figure, we can obviously see the 
impacts of the potential functions V-i,{z) and V-r.{z) on 
the left- and right-handed massive KK modes with even 
parity or odd parity. From the comparison between the 
massive resonant wave functions and massive off-resonant 
ones, we find that the resonances occur only at some par- 
ticular mass eigenvalues, e.g., vn? = 0.072 (even), 0.072 
(odd), etc. for fc = 3, 7? = 0.05, and Vb = 0.05. These 
mass eigenvalues show that they follow the Kaluza-Klein 
parity-chirality decompositions of the massive fermion 
resonances obtained [66]: 

*(x,z) = e-^^Y. [V'L„(x)/(f^z) + V'Rn(^)/r(^)] , 

n 

(48a) 

^{x,z) = [^L„(x)/i°Hz) + 7AR„(x)/fJ(^)] , 

n 

(48b) 

where the superscripts "E" and "O" stand for even parity 
and odd parity, respectively. As we know, this demon- 
strates that Dirac fermions can be composed of the left- 
handed fermions with odd parity and the right-handed 
ones with even parity, and vice versa [64, 66]. Specifi- 
cally, the parities and chiralitics of the resonant massive 
KK modes are conserved for left- and right-handed reso- 
nant massive fermions. 

Besides, we note that the potential of left-handed 
fermions for k = i becomes a double-well potential as 
discussed in Ref. [64, 66], and this can be seen from 
Fig. 5(b). After simple calculations, we find that the 
potential of left-handed fermions is always a double- well 
one for F{(j)) = (j)^ with odd fc > 3, whereas the potential 
of right-handed fermions is a single- well potential all the 
time. However, there exist fermion resonances for both 
kinds of potential wells. 
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FIG. 5. (color online). The shapes of massive resonant KK modes, potential functions, off-resonant wave functions of left- 
handed and right-handed fermions with even parity and odd parity for F{4>) = 4>''^ . The parameters are set to Zmax = 120, 
Vb = 0.05, and rj = 0.05. The terms "res." and "off res." stand for resonant states and off-resonant states, respectively. 
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FIG. 6. The probability Pl,r (as a function of m^) for finding massive KK modes of left- and right-chirality fermions with mass 



around the brane location for F{(f>) 



The solid lines and dashed lines are plotted for the odd-parity and even-parity 



massive fermions, respectively. The parameters are set to Zmax = 120, Vb = 0.05, and rj = 0.05. 



According to the knowledge of quantum mechanics, 
we know that, when a microscopic particle encounters 
a potential well, the microscopic particle will stay in 
the potential well with a limited lifetime. Next, we 
will carefully study the resonant processes. Since the 
Schrodinger equations (37) can be recast into the form of 

cI,rCl,r/l,r(2:) = m2/L,R(2), the probability for find- 
ing the massive KK modes around the branc location 
along extra dimension is |/l,r(2)P- From the previous 
section, we know that the energy density of the thick 
brane concentrates in a finite range along the direction 
of the extra dimension. Strictly speaking, we can inter- 



pret |/l.r(0)P as the probability for finding the massive 
KK modes at the center of the brane. However, for the 
massive KK modes with odd parity, the wave functions 
at the center of the brane are /l'r(O) = 0. It means that 
the probability is poorly defined, and we cannot under- 
stand the resonances very well. Therefore, without loss of 
generality, we define a relative probability to express the 
probability of massive KK modes with even and odd par- 
ities around the brane location. The relative probability 
is defined as [64] 

A ,,,, 
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where the relationship between Zb and ^max is chosen 
as Zmax = lO^b, so the probabihty for the plane wave 
modes with the eigenvalue is 1/10. When wc set the 
same parameters Vq = 0.05 and 77 = 0.05 as in Fig. 5, 
for F{<f)) = (j)'^ with = 3, the probabilities Pl.r (as a 
function of m^) for finding massive KK modes of Icft- 
and right-chirality fcrmions around the branc location 
are depicted in Fig. 6. Wc find that the probability 
is maximized at some mass eigenvalue (a series of res- 
onant peaks) and in the figure the plateau corresponds 
to Zb/^max = 0.1. The number of resonant peaks for 
left-handed KK modes is equal to that of the ones for 
right-handed KK modes. In a sense, this shows that a 
massive Dirac fermion can be composed of the left- and 
right-handed massive KK modes. 

Furthermore, in the numerical calculations, we find 
that the magnitude of the probability for finding a 
fermion around the brane is related to the step size in the 
numerical experiments. Accordingly, in order to analyze 
the fermion resonances, the step size of the coordinate 
used in the numerical calculations must be small enough 
to reflect the true face of resonances. In the subsequent 
discussions, we take a small enough step size to study the 
probability of the resonant states. In this case, more de- 
tailed graphs of resonant probabilities with respect to the 
resonant mass m can be obtained. These resonant peaks 
correspond to the massive resonant KK modes with even 
and odd parity obtained in Fig. 5. We can clearly see 
that more notable details (such as the maximum Pmax of 
resonant probabilities) of resonant peaks are shown com- 
pared with the resonant peaks in Fig. 6. At the same 
time, from Fig. 6, we can see that when the mass eigen- 
value increases, the resonant peaks become thicker and 
thicker. On the other hand, the resonant peaks with 
relatively smaller mass eigenvalues are narrower, and the 
resonant peak with the smallest mass eigenvalue (the first 
peak) is the thinnest one. Later, we will know that the 
narrowest resonant peak has the maximum resonant life- 
time. 



TABLE I. The eigenvalue m^, mass m, widtli F, lifetime r, 
and maximum probability Pmax for resonances of left and 
right chiral fermions with odd-parity and even-parity solu- 
tions for F{(j>) — 0^. C and V stand for chirality and parity, 
respectively. £ and TZ are short for left handed and right 
handed, respectively. The parameters are set to Zmax ~ 120, 
Vo = 0.05, and 77 = 0.05. 

k C V m„ r r 

Odd 0072 0269 0.000 3 3 505.3 0999" 
Even 0.224 0.474 0.003 343.0 0.971 
Odd 0.386 0.621 0.013 78.97 0.731 
Even 0.551 0.742 0.043 22.99 0.362 
Ewn (1072 0269 0.000 3 3 422.2 0998" 
Odd 0.224 0.473 0.003 324.0 0.972 
Even 0.385 0.621 0.014 71.21 0.716 
Odd 0.549 0.741 0.051 19.63 0.350 



To further analyze the lifetime of a fermion resonance, 
first, we define the width F = Am of a resonant state as 
the width at the half maximum of a resonant peak [57] . In 
this case, a massive fermion will disappear into the fifth 
dimension after staying on the brane for some time r ~ 
F^^. Thus, T is called the lifetime of a fermion resonance 
mentioned above. After numerical calculations, we can 
get a lifetime from each peak of the fermion resonance. 
In Table I, we list the eigenvalue m^, mass m, width F, 
lifetime t, and maximum probability Pmax for resonances 
of left and right chiral fermions with odd-parity and even- 
parity solutions for F{<j>) = (j)^ ■ It is interesting that the 
potential wells of the massive chiral fermions are not deep 
enough for fc = 1, Vq = 0.05, and 77 = 0.05, consequently, 
there is no resonant state in this situation. However, if 
we choose a larger 77, the height of the potential of the 
massive chiral fermions will become larger, then a series 
of massive fermion resonant states will appear. 

Note that the data in Table I only refiect the informa- 
tion of the fermion resonances for fc = 3. After further 
calculations, we find that, for a given set of parameters, 
the number of resonant states increases with fc. The res- 
onant mass eigenvalue m of a left-handed fermion with 
odd parity is equal to that of a right-handed fermion with 
even parity. On the other hand, the resonant mass eigen- 
value m of a right-handed fermion with odd parity is also 
equal to that of a left-handed fermion with even parity. 
This means that, in this case, the massive KK modes 
satisfy the equations of the Kaluza-Klein parity-chirality 
decompositions (48). In other words, a massive Dirac 
fermion can be composed of the left- and right-handed 
massive KK modes, as mentioned above. Indeed, this 
result can be understood in the context of supersymmet- 
ric quantum mechanics. The Hamiltonians in (35) and 
(36) can be factorized as Hi, = A and = AA^ . 
Hence and are actually conjugated supcrsymmet- 
ric partner Hamiltonians with the superpartner poten- 
tials Vi,{z) and Vr(z). This leads to the correspondence 
between the spectra of left- and right-handed fermions 
except a normalized zero mode of J-lq^z), i.e. the spectra 
of i?L and are degenerate. The underlying reason 
for the degeneracy of the spectra of Hi, and iJR can be 
understood from the properties of the supersymmetric 
algebra[74]. 

For the lifetimes of the massive resonant states, there 
exists a following rule. For a relatively smaller resonant 
mass eigenvalue m, the resonant state has a relatively 
larger lifetime r. Actually, this rule can be seen obvi- 
ously from the profiles of the peaks of the resonant prob- 
abilities. From Fig. 6, wc can clearly see that, when the 
magnitude of a mass eigenvalue increases, the resonant 
peak would become thicker. The resonant peaks with rel- 
atively huger mass eigenvalues would be broader, and the 
resonant peak with the biggest mass eigenvalue (the last 
peak) is the thickest one. When the resonant peaks be- 
come wide enough, the massive KK modes will not stay 
on the brane, and the time scale r 0. Therefore, the 
massive KK modes will tunnel through the brane at a 
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FIG. 7. The mass spectrum of resonances for F{(f>) = (j)'^. The 
parameters are set to Zmax = 120, Vb = 0.05, and rj = 0.05. 



TABLE II. The eigenvalues m^, mass m, width F, lifetime 
r, and maximum probability Pmax for resonances of left and 
right chiral fermions with odd parity and even parity for 
F(<ji) = 4>. C and V stand for chirality and parity, respec- 
tively. £ and TZ are short for left- handed and right handed, 
respectively. The parameters are set to = 1, Zmax = 30 for 
Vb = 0.1, and Zmax ~ 60 for Vb — 0.5. 



Vb 
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V 


mi 


m„ 
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T 


p 

max 




C 


odd 


0.815 


0.902 


0.004 


266.1 


0.987 


0.1 


even 


1.417 


1.191 


0.047 


21.34 


0.548 


n 


even 


0.814 


0.902 


0.004 


261.6 


0.989 




odd 


1.416 


1.190 


0.049 


20.44 


0.562 


0.5 


c 


odd 


1.547 


1.244 


0.152 


6.569 


0.408 


7^ 


even 


1.522 


1.234 


0.166 


6.021 


0.418 



moment, consequently, they will become free both in the 
bulk and on the brane. On the contrary, for the resonant 
state with the smallest mass eigenvalue, it will have a 
very long time to stay on the brane. The massive KK 
modes with a huge lifetime r will gradually pass through 
the brane, and slowly fade away into the bulk. The ex- 
treme situation is that, when r tends to infinity, i.e., 
r — >■ oo, the resonant state will become a bound state. 

In addition, we can plot the mass spectrum of the 
massive resonant KK modes. For F{(j)) = (j)^ with 
fc = l,3,5,..., the mass spectrum of the massive fermion 
resonances is depicted in Fig. 7 for fc = 3. We can see that 
the number of the resonant states of the massive fermions 
with left- and right-handed chiralities increases with fc. 
However, in the later discussions, we will find that the sit- 
uation is slightly different. The resonant states here arc 
the ones whose resonant mass eigenvalues are below the 
maxima li.,R,max of the potentials (note that VL^max and 
^R,max are the same magnitude). In general, the num- 
ber of the right-handed resonant states of the massive 
KK modes is equal to that of the ones with left-handed 
chirality. 



2. Case II: Vo as a variable 

In this subsection, we mainly discuss the impacts of 
the brane parameter Vq (it decides the thickness of the 
brane by the relation d = 5.074 5/y/Vo) on the resonant 
states of the massive KK fermions, especially on the res- 
onant probabilities, resonant masses, and the number of 
resonant states. We find that the impacts of Vq and k on 
the resonances are different. 

For given 77, fc, and various values of Vq, we study the 
resonances of massive KK modes. As investigated in the 
previous subsection, using the Numerov method, we nu- 
merically solve the Schrodinger equations (37) with the 
potential functions (39). Following the same routine used 
above, for F((/)) = cj), the eigenvalue m^, mass m, width 
r, lifetime r, and maximal probability Pmax are listed in 
Table II for massive resonances of left and right chiral 
fermions with odd and even parities. 



From Table II, we can clearly see that all the values 
of the eigenvalue m^, mass m, width T, lifetime r, and 
maximal probability Pmax f^oUow the same rules obtained 
in the previous subsection. However, when Vq becomes 
larger, the results are quite different from those with 
smaller Vq. In this case, we find that the maximum values 
l^.max and Vn.max of the potential functions of the left- 
and right-handed KK modes are clearly the same mag- 
nitude. The potential for the smaller Vq is deeper than 
that of ones for bigger Vq, while the height of potential re- 
mains unchanged for various values of Vq ■ Consequently, 
as Vq increases, the number of resonant states decreases. 

We find that the probability of the first peak of reso- 
nant states with right-handed chirality is the maximum, 
and it is equal to that of the one with left-handed chi- 
rality. For both chiralities, as the resonant mass eigen- 
value increases, the resonant probability decreases. The 
resonant probabilities of left- and right-handed massive 
KK fermions have a good match at the same resonant 
mass eigenvalue. Accordingly, there is a one-to-one cor- 
respondence between the peaks of left- and right-handed 
fermions at the same resonant mass eigenvalue. As al- 
ready discussed, the mass spectrum can be plotted from 
the data in Table II. Until now, we have seen that there 
exists a good match of the probabilities Pl and Pr for 
all of the excited states and its origin is the consistency 
of the behavior of the potentials and V^. We can see 
that, for each n, which labels the nth resonance, within 
some given numerical error, the resonant mass m„ of the 
resonant KK modes with left-handed chirality is equal to 
that of ones with right-handed chirality. In other words, 
a massive Dirac fermion could also be composed of the 
left- and right-handed massive KK modes. 

In a word, in this section we mainly investigate the 
resonant problem of the massive chiral KK modes with 
odd and even parities. For various values of fc and Vq, 
we obtain the eigenvalue m^, mass m, width F, lifetime 
T, and maximum probability Pmax for the resonances of 
left and right chiral fermions with odd and even parities 
for F{(j)) = (j)^ with odd fc = 1,3. The results can be 
extended to P(c/>) = (t>^ with odd fc = 5, 7, • • • . Besides, 
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we obtained the graphs of the mass spectra of the mas- 
sive KK resonant states. For the two cases, as discussed 
in Refs. [64, 66], the resonant states of the massive KK 
modes satisfy the Kaluza-Klein parity-chirahty decompo- 
sitions. As the parameter Vq becomes larger, the depth 
of the potential wells Vl and Vr becomes shallower, while 
the height of the potentials remains unchanged. Conse- 
quently, as Vo increases, the number of resonant states 
decreases. In other words, when the thickness of the 
brane becomes thicker, the number of resonant states in- 
creases. The resonant probabilities Phimn) and Pr(to„) 
of the resonant massive KK modes have a good match 
at some resonant mass eigenvalue for all of the excited 
states. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have investigated the fermion reso- 
nances in the background of a single scalar field generated 
thick branc with a piecewise warp factor. Using the warp 
factor inspired by the one in Ref. [02], we obtained the 
kinklikc background scalar field without derivative jumps 
under three junction conditions. The distance d between 
the two boundaries, which is just the thickness of brane, 
can be obtained in terms of one free parameter Vq. Next, 
we introduced the interaction between a massless bulk 
fermion and the background real scalar field by means of 
a general Yukawa coupling rj'^ F {(t))'i> in five-dimensional 
space-time, where = 0*"' with odd k = 1, 3, 5, . . .. 

When the parameter k is treated as a variable, we 
found that the coupling between the scalar and fermion 
infiuences the resonances. The number of resonant states 
of the massive KK fermions with the left- and right- 
handed chiralities increases with k. Here, the number 
of resonant states is defined as the number of those res- 
onant states with eigenvalue lower than the highest 
point of the potential. We found that the number of the 
resonant KK modes with left-handed chirality is equal 
to that of the ones with right-handed chirality. This 
means that the fermion resonant states of massive KK 
modes satisfy the Kaluza-Klein parity-chirality decom- 
positions [GG]. For every k = {1,3,5,...}, the resonant 
KK fermions with smaller mass eigenvalues have larger 
lifetimes than those with bigger mass eigenvalues. On 
the other hand, we can see that the wider resonant peaks 
have smaller lifetimes. Furthermore, the first resonant 
peak is the narrowest one, and it has the longest life ex- 
pectancy. Accordingly, the KK fermion resonances with 
smaller mass eigenvalues can stay on the thick brane for 
a longer time. In other words, these quasibound states 
(resonant states) are more stable. Conversely, the KK 
fermion resonances with huger mass eigenvalues can only 
stay on the brane for a little time, i.e., the resonances 
become too unstable and cease to appear [G3]. 

When the brane parameter Vq is treated as a variable, 
we found that the thickness of the brane also infiuences 
the fermion resonances. The number of the peaks of res- 



onant probability of the massive KK fermions with left- 
and right-handed chiralities decreases with Vq, namely, 
the number of resonant states decreases with Vq. When 
Vo varies, the maximum values ^L.max and VR,max of 
the potential functions of the left- and right-handed KK 
modes are the same magnitude. Accordingly, the proba- 
bilities of the resonant KK modes of left-handed fermions 
arc equal to that of the ones of right-handed fermions. 
Further, the masses and lifetimes of left and right chiral 
resonances are almost the same, which demonstrates that 
it is possible to compose a massive Dirac fermion from 
the left and right chiral resonances. In this case, the res- 
onant probabilities of the left- and right-handed resonant 
KK modes with the same mass are the same for all the 
massive resonant KK modes. There is a good match for 
the resonant probabilities of the fermion resonances with 
left and right chiralities. 

Indeed, if a resonance is supposed to represent a 
metastable four-dimensional massive fermion field, it 
must contain both left and right parts, which are mixed 
by the four-dimensional Dirac equations [see Eqs. (28) 
and (29)]. We can see that there is no mismatch for 
peaks of the probabilities Pl and Pr (as a function of m^) 
for finding massive KK modes of left- and right-chirality 
fermions. The numbers of peaks of Pl and Pr are equal 
to each other. So the resonant masses of the left- and 
right-handed fermions must be equal to each other. So 
our computations are consistent with the coupled four- 
dimensional Dirac equations. Actually, the definition of 
the number of resonant states in Ref [GG] can be used 
because there is a good match for the maximum values 
^L,max and Vn^niax of the potential functions of the left- 
and right-handed KK modes. From the definition, we 
can see that only the resonant states with eigenvalue 
lower than the highest point of the potential are consid- 
ered as resonant states because they have obvious reso- 
nant peaks, namely, they have relatively larger resonant 
probabilities. In other words, the resonant states with 
unobvious resonant peaks were neglected. So the reso- 
nant states with larger mass eigenvalues (the ones with 
quite small resonant probabilities) were not taken into 
account. The prior condition is that the maxima of the 
potential functions Vl and Vr are equal to each other, 
or at least they have the same magnitude. Actually, the 
numbers of resonant states can be considered as numbers 
of peaks of Pl or Pr . Although the resonant probabilities 
(whose resonant masses is larger than the threshold value 
of Va^niax) are quite small, we should take them into ac- 
count. Thus, from the physical point of view, all the res- 
onant peaks describe the resonant states. Essentially, the 
resonant behavior coincides with the qualitative change 
of the energy density of the brane. 

With respect to resonances for chiral fermions on the 
thick brane with a piecewise warp factor, we can resume 
our conclusions with the following points: (i) The thick 
brane with a piecewise warp factor was derived under 
three junction conditions which ensure the correct Israel 
junction conditions. The kinklike scalar field without 
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derivative jumps can be given. Especially, the continuous 
scalar potential and energy density were also obtained, 
(ii) The thick brane without internal structure also fa- 
vors the appearance of resonant states for both left- and 
right-handed fermions. When the brane becomes thicker 
(Vo becomes smaller), the number of resonant states in- 
creases, and vice versa, (iii) The scalar-fcrmion coupling 
influences the resonant behaviors of the KK fermions. 
The resonant properties (such as resonant probability, 
lifetime, mass, and number) of the left-handed fermions 
are the same as the ones of the corresponding right- 
handed fermions. (iv) Strictly speaking, all of the res- 
onant peaks describe the resonant states. The definition 
of the number of resonances given in Ref. [66] should be 
redefined. Indeed, the number of resonant states can be 
considered as number of resonant peaks. The fermion 
resonant states of the massive KK modes satisfy the 
Kaluza-Klein parity-chirality decompositions, (v) It is 
worth noting that many systems with boundaries cannot 
be investigated in the thick brane scenario, instead the 
thin branes as the boundary terms should be added into 
the action. Such systems can be studied with the method 



in Ref. [19]. The properties of the fermion resonances in 
thick brane models and the problems in the brane with 
boundaries are rich and interesting. We leave these issues 
for future works. 
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